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(57) Abstract 

The present invention of interfering continuous-wave signals are automically detected in the demodulated output of a quadrature 
amplitude modulation data communications receiver. The transmitter (16) is re-tuned in ordw to avoid the detected mterferer. Statistiral 
information concerning the received signal is accumulated at a receiver (18) and is analyzed to determme whether a reducUon m the 
received bit-enir-raie is due to the presence of an interferer. If an inteferer is determined to be present widun the signal transmission 
bandwidth then the potentially available transmission band is spectrum analyzed to find a new portion of that band m which the signal can 
be transmitted free of tiic interferer. The results of die statistical analysis and die spectral analysis are conveyed back to the transmitter | 
(16) which then readjusts the transmission parameters and forwbds the new parameters to the receiver (18). The system resumes operation 
in the new transmission band with no interference and with die tit-error-rate restored to its required level. 
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INTERFERENCE DETECTION AND AVOIDANCE TECHNIQUE 

Technical Field 

The present invention relates generally to data 
5 communications, and more particularly, to a method and 
apparatus for automatically detecting the presence of an 
interfering signal in a data communications receiver and 
retuning a transmitter to avoid the detected interfering 
signal . 

10 

Background Art 

To meet the growing need for ever increasing digital 
data bandwidth for new subscriber services (e.g., high- 
data-rate internet service, video telephony, and high 
15 definition television) telephone companies are looking 
towards the use of very-high- speed digital subscriber 
lines (VDSL) . VDSL provides a means to carry such data 
into the home or business over the existing copper wires 
used by plain old telephone service (POTS) . VDSL 
20 carries digital data at bit rates reaching 52 Mbps using 
carrierless AM/PM (CAP) modulation, which is related to 
quadrature amplitude modulation (QAM) . Modulation 
contellations can reach 256 symbol points with symbol 
rates up to 6480 Kbaud. The CAP modulation 
25 constellations look the same as the corresponding 

baseband QAM constellations. The signals are carried 
within the spectral band up to 30 MHz, which is above 
the frequency band used by POTS on the same wire pair. 
Due to the high loss of the copper telephone wires at 
30 these frequencies, the VDSL signals are carried on the 
telephone wires only over the "last mile" (i.e., the 
last segment of copper wire between a central office anc 
the user's premises) . ^ 
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Within the spectral band of 0.5 MHz to 30 MHz, the 
VDSL transmissions over the telephone wires are subject 
to radio frequency interference (RFI) from AM radio 
broadcasts and ajmateur radio transmissions, which can be 
5 picked up by the telephone wires acting as receiving 

antenna. This interference can corrupt the received and 
demodulated CAP/QAM signals, causing symbol and bit 
errors in the received data. The AM broadcast signals 
remain indefinitely at the same frequencies in any 
10 neighborhood within range of the local AM radio 

transmitters. The amateur radio signals, however, will 
vary in frequency throughout the day as the operator 
tunes his transmitter to optimize the HF-band skywave 
propagation conditions over the communication links he 
15 is attempting to maintain. There are several narrow 

spectral bands allocated to amateur radio within the 0.5 
MHz to 30 MHz band carried by VDSL transmission lines. 
A VDSL transmission line will normally be interfered 
with only by amateur radio transmitters within a few 
20 hundred feet of the transmission line. 

Narrowband interferers, such as AM radio broadcasts 
and amateur radio signals, are not the only inter ferers 
that corrupt the CAP/QAM signals. Wideband noise that 
extends across a large portion of the VDSL spectrum may 
25 also interfere, causing undesirable changes in the 
carrier-to-noise ratio (CNR) and the bit error rate 
(BER) . Wideband noise can be caused, for example, by 
electrical machinery, internal combustion engines (e.g., 
lawn mowers) or fluorescent lights operating in the 
30 immediate vicinity of a cable drop to the premises. 

Narrowband interferers tend to be more predictable and 
hence somewhat easier to avoid, for example, by 
detecting the frequency of the interferer and subsequent 
adjustment of the CAP/QAm VDSL signal for avoidance 
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thereof. Wideband interferers, on the other hand, are 
somewhat more random, and cannot generally be accounted 
for prior to data transmission and are therefore more 
difficult to avoid. Especially difficult is detection 
5 and characterization of the interferers without 
interruption of the data signal for analysis. 

Thus, what is needed is a method and apparatus for 
recognizing the difference between narrowband 
interferers and wideband noise sources without relying 
10 upon spectrum analysis of* the" received signal except to 
locate a new transmission band. This is accomplished by 
the present invention whereby a determination is made of 
the presence or absence of a narrowband or continuous 
wave (CW) interferer, and in many cases inferring its 
15 carrier frequency without the need for interrupting the 
transmitted data stream. The transmitted data stream is 
only interrupted when an inteferer is detected and a 
spectrum analysis is needed to find a new transmission 
band. The interferer presence determination can 
20 generally be accomplished when the interferer level is 
low enough that the presence of the interferer cannot be 
inferred by visual examination of the constellation plot 
of the demodulated signal. When a narrowband signal is 
present, the shape of the symbol clusters in the 
25 constellation plot changes from a centrally peaked 
distribution, characteristic of a signal in Gaussian 
noise, to a ring as shown in FIGs . lA and IB. The 
present invention recognizes this transition 
automatically and at low interference levels where the 
30 transistion to a ring is not yet visually apparent. 

Adjustment of the transmitted symbol rate, bandwidth and 
carrier frequency is provided in order to avoid the 
troublesome interferer and other interferers that may be 
present in the signal environment. 
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Disclosure of Invention 

According to a first embodiment of - the ipresent . 
invention, a method of non-intrusively detecting whether 

5 an interference signal is present with a received signal 
of interest is provided. The interference signal is 
detected by accumulating statistical information 
associated with the signal of interest. - The statistical 
information is analyzed to determine whether an 

10 interference signal is present and is intrusive. If an 
intrusive interference signal is present, .it is analyzed 
to determine its frequency. An alternate*^ spectral 
region is next identified that is free of interferers. 
Lastly, communications are adjusted by re-tuning or 

15 modifying data rates to allow continued communications 
without intrusive interference. 

In another embodiment of the present invention, a 
modem capable of non-intrusively detecting an unwanted 
interferer while a receiving of a data signal is 

20 provided. The modem includes a transmitter for 
transmitting information to a second modem and a 
receiver for receiving the data signal from the second 
modem. A memory is connected to the receiver for 
storing statistical information associated with said 

25 data signal, including data sample values in each 

received symbol cluster in the data signal's amplitude 
phase constellation. A processor measures constellation 
statistics from the statistical information, determining 
a carrier-to-noise ratio and bit error rate of said data 

30 signal, determining whether the bit error rate exceeds 
an acceptable limit by testing whether symbol cluster 
changes from a Rayleigh distribution to a Ricean 
distribution, and when the bit error rate exceeds said 
acceptable limit, interiipting the received data signal 
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to spectrum analyze a transmission band to find an 
interference free region of the transmission band and 
communicating the interference free region to the second 
modem for adjustment of transmission parameters. 

5 

Brief Description of Drawings 

FIGs. lA and IB are an illustration of the effect of a 

sinusoidal tone interferer upon the symbol clusters 
in a quadrature amplitude modulation (QAM) - 
10 constellation. 

FIG. 2 is a block diagram of a communication system 

using the present invention. 
FIG. 3 is a flow chart of the processing step of 
measuring the constellation statistics. 
15 FIG. 4 is a flow chart of the processing step of 

determining the carrier-to-noise ratio (CNR) and * 
bit error rate (BER) . 
FIG. 5 is a graph of the standard deviation of 

determined CNR versus size of the input data set 
20 for a 16-QAM signal with a CNR of 12.8 dB. 

FIG. 6 is a graph of the standard deviation of 

determined CNR versus size of the input data set 
for a 64-QAM signal with CNR of 19.1 dB. 
FIG. 7 is a graph of the ratio of the mean to standard 
25 deviation, pi(z), versus z for a Rice distribution. 

FIG. 8 is a graph of the skewness, p3(z), versus z for a 

Rice distribution. 
FIG. 9 is a graph of the kurtosis, p4(z), versus z for a 
Rice distribution. 
30 FIG. 10 is a graph depicting cluster skewness as a 

function of CNR and interferer amplitude versus 
inteferer power. 
FIG. 11 is a graph depicting cluster skewness as a 

function of CNR an^d interferer amplitude versus 
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normalized interferer amplitude . 

Bes^ Mode for Carrying Out the Invention ^ : : 

The present invention provides- an improved method 
5 for recognizing the presence of a CW interferer in the 
demodulated output of a CAP or QAM digital data 
receiver. The method relies on the fact that the vector^ 
combination of the interferer with a- CAP /QAM symbol' - '^ 
cluster changes the radial distribution of sample 
10 locations in the cluster from a Rayleigh distribution in 
the presence of only Gaussian noise to a Ricean" 
distribution in the presence of such an interferer. The 
invention uses statistical analysis of the positions of 
a set of samples in the demodulated symbol constellation 
15 to decide whether the distribution is sufficiently 

shifted towards Ricean to declare that an interferer is 
present. 

Three different measures of the shape of a 
probability distribution function were tested: the ratio 

20 of the mean to standard deviation, the skewness or third 
central moment, and the kurtosis or fourth central 
moment. Of these measures, the skewness was found to be 
best correlated with the received carrier-to- 
interference ratio (CIR) , and the kurtosis was also 

25 found to provide a useful degree of correlation. By 
calculating either the average skewness or average 
kurtosis over all clusters in the constellation, and 
testing this average skewness or average kurtosis 
against a threshold, the present invention determines 

30 that an interferer is present if the average skewness or 
average kurtosis is below the threshold, and that no 
interferer is present if the average skewness or average 
kurtosis is above the threshold. 
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While a fixed threshold is possible, it has been 
found that at high CNRs the interference levels 
(measured as a fraction of the noise margin), that can be 
positively recognized in this way are higher than the 
5 interference levels that can be recognized at low CNRs. 
In high noise (low CNR) situations there is not much 
difference in the, skewness or kurtosis of an 
interference-free signal and one that has a moderate 
level interf erer.^that could otherwise be detected in a 
10 high CNR environment.. With a small set of sample data, 
the statistical fluctuations in the computed skewness or 
kurtosis can be , larger, than the significant differences 
in skewness. However, with a large enough sample set-, 
even small differences in skewness or kurtosis and con- 
15 sequent CIR can be detected, even at low CNR, providing 
that the threshold criterion is made a function of the 
CNR (which is routinely measured) . Generally, skewness 
has been found to be more stable with respect to 
sampling fluctuations than is kurtosis. The main 
20 limitations will be the speed with which the data 
samples can be transferred to the data processing 
computer and the processing speed of the computer. For 
a cost-limited unit with limited computing resources, 
such as a VDSL subscriber modem, there may be design 
25 limitations on the number of data samples or the 
sampling rate that can be handled. 

According to the present invention, the receiver 
accumulates certain statistics regarding the data sample 
values in each received symbol cluster in the signal's 
30 amplitude -phase constellation. These statistics are 
used to calculate the CNR (actually, the carrier to 
noise-plus- interference ratio) and the BER for the 
received signal constellation. If the BER falls below 
an acceptable limit, then the average cluster skewness 
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and/or kurtosis is also computed for the constellation. 
A table 'lookup provides the threshold- skewness and 
kurtosis- values at -the calculated CNR below which a 
calculated skewness- or t. kurtosis would indicate the 

5 presence of an . interf erer of amplitude sufficient to 
significantly affect the measured CNR. The calculated 
skewness and/or kurtosis are then compared to the 
selected threshold skewness and/or^ kurtosis , and if the 
calculated skewness and/or kurtosis are below the 

10 threshold;'- then the presence ofi-an interf erer is 
declared. - :^ ' , - 

The frequency of the interferes is then determined 
by spectral analysis of the received transmission band 
for the signal of interest using, most commonly, an FFT. 

15 If the interferer is strong (small CIR) , the spectral 
line of the interferer may show up above the spectrum of 
the SOI. However, if the amplitude of the interferer is 
approximately the same as that of the SOI, it will 
generally be necessary to interrupt transmission of the 

20 SOI long enough to collect interference data to spectrum 
analyze. In some cases, it is possible to infer the 
frequency of the interferer if there can only be one 
such frequency. For example, in VDSL transmission there 
may be only one allocated amateur radio band within the 

25 spectral band used for a particular VDSL transmission in 
which interference is detected. This will then be the 
interferer frequency. It may then be possible and 
acceptable to avoid the interferer by reducing the 
transmission bandwidth .{and the corresponding data rate) 

30 to the segment of the original transmission band on one 
side of the interferer, so that the interferer is then 
outside the new VDSL transmission passband. 
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If reducing the data rate is not an acceptable 
alternative, _ then-, it^ will be necessary to find, another 
passband within -the -available VDSL., spectral range that 
is free, of - interferers..,. This requires spectrum^ ...^ 

5 analyzing the entire VDSL spectral range, not just that 
passband in current use. This must be done before any 
bandpass filtering. in the- receiver, and will again , 
require interrupting;,, the , VDSL data transmission long 
enough to collect data . for spectral analysis of. ^the 

10 interference environment;- The advantage of theVpr.esent 
invention is that this interruption of the SOI data 
transmission -need only be done when it has been -^- 
determined from constellation analysis that there is an 
interf erer present which is causing unacceptable 

15 degradation of the BER. If the BER is so degraded, then 
it does not matter that data transmission is 
interrupted, since the received data already has too 
many errors. It is better to interrupt transmission and 
find an acceptable transmission band so that suitably 

20 error- free transmission can be restored as quickly as 
possible. 

The present invention detects when an interferer is 
present that is degrading the BER of the SOI. The 
invention then provides means for determining a 

25 transmission band that will be free of degrading 

interferers and then retuning the transmitter to the new 
transmission band. 

Referring now to the figures, and for the present 
to Fig. 2, therein is shown a block diagram of the 

30 present invention. In the transmitter 16, the input bit 
stream 20 enters the symbol encoder 22 where the input 
bit stream 20 is encoded into QAM/CAP modulation. The 
resulting symbol sequence has m bits per symbol in a 
constellation of M symbols. Generally, m=log:,M. In the 
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modulator 24, the symbol sequence out of the symbol 
encoder 22 is modulated onto' in-phase and quadrature - r 
rails at ^a"- carrier- frequency derived ^from .carrier, 
oscillator 26. The' modulated signal-is then boosted in - 
5 power by the power amplifier 28 and the signal of 
interest (SOI) is transmitted over the forward 
transmission channel 30 to the receiver 18. Along the 
forward transmission- channel 30,- external radio- 
frequency interf erence-32 is added to the transmitted 
10 SOI, the addition being modeled by the summation 34. 

The combined signal of ^" SOI plus RFI enters "the 
receiver 18 where' it is first demodulated in the 
demodulator 36. This includes locking to the carrier 
frequency and symbol clock rate, separating the in-phase 
15 and quadrature components, downconverting from the 

. carrier frequency to baseband, and sampling the signal 
at the symbol clock rate. The output of the demodulator 
36 is a sequence of in-phase and quadrature component 
values for the sum of SOI plus noise plus interference 
20 at each symbol sample. The output of the demodulator 3 6 
is passed through the equalizer 38, which compensates 
for and removes the effects of intersymbol interference 
and changes in sample values caused by channel 
distortion, filtering, and mult ipath propagation (e.g., 
25 extraneous paths involving bridged taps, i.e., open 
circuited twisted pairs connected in shunt with the 
main-path twisted pair on a VDSL line) . The equalized 
signal goes to slicer 40 which tests the in-phase and 
quadrature values against decision boundaries and makes 
30 a decision as to the symbol state for each equalized . 
sample. The symbol stream 42 out of the slicer then 
goes to the symbol decoder 44. The symbol decoder 44 
converts each symbol to the binary bit sequence 
corresponding to the input symbol state and outputs a 
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binary bit stream 46 representing the information 
conveyed by^the- * system The . bit-, ,s tream^ 45 goes^ to . an _ ..^ 
interpreter (not shown) .which separates. the payload _ 
information from the overhead control bits. . . 

5 The signal at the output of the equalizer 38 is 

sampled by the left hand pole 66 of- double-pole .switch 
64, and the signal- at the output of^ the- slicer 40 is. 
sampled by^^the right hand pole 68^of double-pole switch 
64. Both poles of-, the- switch 64 are sampled in^ _ - : 

10 synchronism under-_.control . of .receiver, controller 62 . . 
Both sets of - sample, data are.^stored^in data storage 
buffer 72. within computer 76, with the output ..samples of 
equalizer 3 8 being delayed so that -they -are aligned with 
the corresponding output samples of the slicer 40 for 

15 the same symbols in parallel columns of the buffer 72, 
These sample times are determined by control bits in the 
output bit stream 46 and by constraints imposed by the 
data collection capabilities and processing rates of the 
signal processing path within computer 76. For example, 

20 the switch 64 might be commanded to sample only every n- 
th symbol out of the equalizer 38 and its corresponding 
sliced value out of slicer 40. In addition, as will be 
described subsequently, switch 64 will normally be 
commanded to collect blocks of sample data only at 

25 specific times as will be described. 

The data storage buffer 72 consists of two blocks 
of random access memory (RAM) within computer 76. The 
two blocks of RAM are "pingponged" , that is, while data 
is being loaded into one block, the data in the other 

30 block is being processed. When the data in the one 

block has been processed, the contents of that block are 
overwritten with new data and the data in the other 
block are then processed. As is known by those skilled 
in the art, other types* of storage can be used in place 
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of the RAM, including, for example, EPROM, SDRAM, hard 
disk, etc. : • • ^-^ - 

At system startup, dataware, initially transferred.,:, 
over the forward transmission channel 30 using a.^ 
5 relatively low order constellation, say 16-QAM, with a 
relatively low symbol rate and with a bandwidth and . . 
carrier frequency selected to have a high probability of 
being free of ^RFI. ^At startup,, the receiver 18- is- 
automatically set by the receiver, controller 62 to .tune 
10 to the standard startup carrier - frequency , symbol rate, 
and constellation size. Thereafter, prior to any time, 
that the 'transmitter 16 plans- to change either the .-., 
carrier frequency, symbol rate or constellation size, 
the transmitter 16 sends a transmission change 
15 notification to the receiver 18. The receiver controller 
62 extracts the new transmission parameters from the 
output bit stream 46 and sends an acknowledgment over 
reverse control channel 58 to the transmitter controller 
60. Only when the transmitter controller 60 has 
20 received this acknowledgment does it command the 

transmitter 16 to change the transmission parameters. 
At this point the receiver controller 62 has commanded 
the receiver 18 and the computer 76 to expect the new 
transmission parameters. Knowledge of the new carrier 
25 frequency and symbol rate permits the demodulator 3 6 to 
rapidly acquire lock to these new parameters, while the 
slicer 40 decision boundaries are reset to correspond to 
the new constellation. 

Once initial receiver synchronization has been 
30 received and acknowledged to the transmitter controller 
60 over the reverse control channel 58, the transmitter 
controller 60 informs the receiver controller 62 over 
the forward transmission channel 30 and then turns off 
power amplifier 28 to interrupt the transmitted signal: 
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The receiver carrier frequency and symbol rate phase 
locked. loops, in demodulator 36 continue to flywheel at 
the previously- acquired, carrier frequency and symbol 
rate during this, interruption. - Also during the 
5 interruption, switch 70 is closed and switch 64 is 
opened, so that the radio frequency, interference 
environment, is sampled and digitized by analog-to- ^ ^ - _ 
digital; converter 78._ Computer 76 then collects 
digitized time samples of the interference environment 
10 which it spectrum analyzes using an FFT 74. The , - 
resulting interference- environment spectrum -is -trans- ^ 
mitted over reverse control channel -58 . to transmitter 
controller 60, which selects a new carrier frequency, 
symbol rate and constellation size that are adequate to 
15 carry the desired pay load data rate while avoiding any 
RFI that may be present in the environment. After 
notifying the receiver controller 62 of the new 
parameters and receiving an acknowledgment, the 
transmitter controller 60 resets the parameters in the 
20 symbol encoder 22 and carrier oscillator 26, and turns 
back on the power amplifier 28 to begin data 
transmission. 

The first processing step on any block of input 
data out of data storage buffer 72 is to measure the 
25 constellation statistics 48. The flow chart for this 
processing step is shown in Figure 3. The computer 76 
receives information in step 302 from the receiver 
controller 62 as to the size of the symbol 
constellation. A block of data storage is then 
30 allocated to each symbol cell at step 304, and the 

symbol data out of the equalizer 38 are sorted into the 
appropriate symbol cell locations at step 306 according 
to the corresponding cell identifiers that come from the 
output of the slicer 40. After the presliced symbol 
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data is sorted into bins by cell, the number of symbols 
in each cell are: counted (step .308), and the.isample mean - 
values and sample, variances of- the. data in . the_„in-phase . 
(x) and quadrature (y) . directions are calculated . for . 

5 each bin at steps 310 and 312. This may be done for 
each symbol cell -in the constellation or, for large 
constellations such as 256-QAM, for just the cells in 
one quadrant. With M cells and sample symbols in*' 
each cell; the sample" means and sample variances " are - 

10 given by: - - - ~ " - ' - 

(1) Sample mean x-position of sample set in cell m: 

(2) Sample mean y-position of sample set in cell m: 

ym~Tj~Svmn' 

15 (3) Sample variance on x-positions of sample set in 

cell m: 

1 



2 



(4) Sample variance on y-positions of sample set in 
cell m: 

1 r^. 1 

20 ^.m=N^^ 

The numbers of samples Nm i^i each cell are counted 
at step 3082, since these numbers are needed in the 
above calculations. 

The next processing step (see block 50 of FIG. 2) 
25 is the calculation of the received signal carrier-to- 
noise ratio (CNR) and probability of error, or bit error 
rate (BER) . The CNR is determined using the method of 
symbol cluster size an'^lysis, wherein the CNR is 
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determined as a function of the constellation closure 
ajd, where is the RMS noise voltage (equal to the RMS 
symbol cluster size) and d is the'noise margin, i;e., 
the minimum (perpendicular) distance," in volts, from the 

5 ideal symbol point to the nearest decision boundary. The 
quantity is called the cluster variance. This method 
is usable for high enough CNRs that each symbol cluster 
is entirely within its decision cell. A flow chart for 
the procedure is shown in Figure 4. First at step 402, 

10 the RMS cluster size is calculated using the 
following equations: 

M 

varx= X(Nm-l)Sx.m. 



15 



1 ^ 



Ivarx + 
rn=y J 



vary 



Then at step 404 the reciprocal of the constellation 
closure d/an is calculated. The CNR is then obtained at 
step 406 as 



CNR(dB)=201ogio 

20 where L is given for each constellation modulation in 
the following table: 
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Modulation 


L 


a 


b 


QPSK - 


2 


2 


1/2- 


.16-QAM 


10 


3 ' 


3/4 


32-Cross QAM 


20 


13/4 


11/13 


64 -QAM 


42 


7/2 


7/8 


128-Cross QAM 


82 


29/8 


53/58 


256 QAM. 


170 


15/4 


15/16. 


9-QPR- . .. 


5.3333-_ 


8/3 


2/3 


25-QPR - . 


16 : 


16/5 


4/5 


49-QPR--" - : 


32 


24/7 


6/-7 __ 


SlrQPR .. . 


53.3333 


32/9 


8/9 



TABLE 1 



The symbol error rate (SER) is obtained at step 408 
5 by first calculating the quantity Q: 

where erfc denotes the complementary error function. A 
lookup table for ef rc ( ) can be stored in the memory of 
computer 76 and interpolated for any value of the 

10 argument. The SER is then calculated using the following 
equation with the appropriate coefficients, as listed in 
the above table: 

SER = agi-bQ. 

As an alternative to using a lookup table for 

15 erfc, the function Q{u) can be evaluated by the 
polynomial approximation 

where 

l + pu 
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u=cJ(J^, and p and Ci,,..,C5 are constants given by 
p = 0.2316419 " 
Cl = 0.319381530'' 
C2 = -0.356563782 - ' 
5 C3 = 1.781477937 

C4 = -1.821255978 ' . ~' " 

Cs = 1.330274429 -/^ ' '~ _- 

The inherent accuracy:, of this approximation is better 
than 1% down to an error probabil:ity-of 10'^^ and 10% 
10 down to - * ~ 

io-^4«. - " ' ' 

Next, the BER is calculated -at-step 410. For full 
response QAM and cross-QAM constellations having M 
symbol cells, the BER is related to the SER by 
15 BER = SER /bits per symbol, 

where 

bits per symbol = log2M. 

To use the constellation closure to determine CNR, 
and from it the SER and BER, it is necessary first to be 

20 able to resolve the individual symbol clusters, and 
second, to be able to determine the cluster variance 
with reasonable accuracy. It is not possible to resolve 
the individual clusters when the RMS cluster size CT^ is 
close to the noise margin d, that is when ajd^l. In 

25 practice, it is not possible to resolve the individual 
symbol clusters when ajyd>K, where K is some number in 
the range 0.5 < K < 1.0. When aJd^K, then 

(CNR),B = 101og,o(iy2K2). 

This gives the limiting minimum value of CNR below which 
30 it is not possible to meaningfully determine the CNR for 
L corresponding to a particular constellation size and 
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for an assumed value of K. These limiting values of CNR 
are listed in the following table 2: 



Modulation 


L 


CNR(dB) 


K = 0.25 


K = 0.5 


K = 0.75 


K= 1.0 


QPSK 


2 


12.0 


6.0 


2.5 


0.0 


16-QAM 


10 


19.0 


13.0 


9.5 


7.0 


32-Cross QAM 


20"' 


22.0 ""^'"^ 


16.0 


12.5 " ' ■ 


ib.o" . 


64-QAM 


42 


25.2 


19.2 


15.7^^ 


13.2 


128-Cross QAM 


82 


28.2 


22.1 


18.6 


16.1 


256 QAM 


170 


31-3 


25.3. 


21.8 


19.3 , , 


9-QPR 


5.3333 


16.3 


10.7 


6.8 


4.3 


25-QPR 


16 


21.0 


15.0 


11.5 


9.0 


49-QPR 


32 


24.0 


18.0 


14.5 


12.0 


81-QPR 


53.3333 


26.3 


20.3 


16.8 


14.3 



5 TABLE 2 

To achieve reasonable accuracy in determining the 
CNR from cluster variance analysis, a large enough 
sample symbol set must be used to reduce the expected 

10 uncertainty, or standard deviation, of the calculated 
CNR to a reasonably small value. Specifying or 
measuring the CNR to an accuracy of much less than 0.1 
dB is generally not meaningful. Therefore, a large 
enough sample set must be used to assure that the 

15 uncertainty in the calculated CNR is no larger than 0.1 
dB. Some Monte Carlo simulations have been done with 
data sets of various sizes to estimate the standard 
deviation of the calculated CNR for various size 
constellations for a CNR that is theoretically adequate 
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to achieve a BER of 1x10"^ (the maximum allowed error, 
rate for a VSDL"modem)~ " ^ - ' 

FIG. 5 shows the standard deviation of calculated 
CNR versus size of the input "data set for a 16-QAM 

5 signal with a CNR of 12.8 dB' (needed to achieve a BER' of ^ 
1x10"^), and Figure 6 "shows the "standard deviation of 
calculated CNR versus, .size of the input data set for a 
64-QAM signal with a: CNR of,19,i dB (needed to achieve a 
BER of 1x10"^). In each case , ^ the standard deviation was 

10 calculated using a set. of ten simulations at each input 
data set size, . which accounts for the still jagged 
nature of each graph. Figure 5 shows that for 16-QAM at 
least 1500 symbol samples are needed, and for 64-QAM at 
least 3000 symbol samples are needed to achieve a 

15 standard deviation of 0 . 1 dB in CNR at the critical CNR 
level. A similar analysis was done for 256-QAM, which 
needs a CNR of 25.15 dB to achieve a BER of 1x10"'^; it 
was found that 50,000 samples were needed to achieve a 
standard deviation of 0.4 dB in the . calculated CNR, and 

20 it is estimated that 0.2-0.3 million samples would be 
needed to achieve a BER of 1x10"^. Consequently, there 
will be a longer time needed between updates of the 
estimates of CNR and BER for 256-QAM than for smaller 
constellations . 

25 For the step of testing the BER (52 in FIG. 2), the 

BER that is calculated using the above described 
procedure is compared to the maximum allowed threshold 
BER. If the calculated BER is at or below the threshold, 
no change is made to the signal transmission parameters, 

30 . and a new set of data is accumulated to update the 

estimated BER and repeat the test. If the calculated 
BER is above threshold, this suggests that there may be 
an interference source '-that is raising the error rate. 
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To check this, the computer 76 goes to the step of 
calculating the average- cluster skewness and/or kurtosis 
54, rather than updating the estimated BER.-rit: must be.^ _ _ 
noted that the sample data in either RAM block of -^the 
5 data storage buffer 72 must not be overwritten with new 
data until . the testing of the BER ,52 is completed, since 
the data, stored in that -RAM block will, be .needed , in the 
subsequent calculation .of skewness and/or_kurtosis 54 if 
the HER test 52 indicates that the BER has risen to an^. 
10 unacceptably high value. To assure this, whenever the 
BER test 52 registers above threshold, the computer 76 
sends a signal informing of this condition to the 
receiver controller 62. The receiver controller 62 then 
opens both poles of switch 64 to prevent new data from 
15 entering into data storage buffer 72. 

When there is no interference present, the x and y 
components of samples in cluster m are each Gaussian 
distributed due to noise. If there are no distortions, 
X and y correspond to the I and Q components. 
20 Furthermore, if there is no gain mismatch, the x and y 
distributions will have the same variance. In that 
case, the radial distance to the samples from the center 
of the cluster will have a Rayleigh distribution. 

When there is an interferer present (assumed to be 
25 a sinusoidal tone at other than the digital carrier 
frequency) , the sinusoidal interferer causes each 
cluster to be ring shaped due to the addition of a 
randomly rotated interference vector of constant 
magnitude to each sample, as shown in Figure 1. If 
30 there were no noise on the system, the interferer would 
transform each ideal constellation point to a sharp 
circle of zero thickness. On the other hand, if the 
interferer had finite amplitude modulation, the locus of 
the interference vector would have a finite thickness. 
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This thickness would add to the width of the noisy, but 
interference-free cluster to result -in an. interference- - 
plus-noisy-symbol cluster-- ring thicker than the . .ring 
that would result from a- pure tonal interferer. - , r-- 

5 When an interferer is present, the radial distance 

to the samples from the center of the cluster will have 
a Rice distribution. - Thus,, to .automatically determine • ' 
whether or not an- interferer is present, and affecting - - 
the digital modem performance, it . is necessary, to 

10 automatically distinguish: between a Rayleigh and a. Rice 
distribution, with, the" same peak. radial distance on the 
Rice distributions for all symbol clusters. For the ._ 
Rayleigh distribution, the probability density function 
is given by 

R.v.c^r)=-^e'^'-. r>0. 
15 ^ 
where 

r = V(Xmk-5<m)^ + (ymk- 7m)^ 
is the radial distance from the centroid O^^Ypn) of 
cluster m to the k-th sample in the cluster, located , at 
20 (Xmkymk)- Also the standard deviations of the cluster 
distribution in each coordinate are assumed equal, so 
that a = ff, =CTw . Using the same notation, the Rice 
distribution is given by 

p^M^-^e' I ^-^J, r>0. 

25 where Io(x) is the zero-th order modified Bessel function 
of the first kind, and a the amplitude of the 
interfering signal. Here, a is the amplitude of the 
interferer. Note that when a = 0, lo(0)=l, so that 
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To decide which type of distribution is present, 
and thus whether or not an interferer is present, 
various moments and functions of moments may, be 
calculated for the sample distribution for each cluster, 
5 and these moment functions may be compared .with the - 
theoretical values for the Rayleigh and Rice 

distributions. This will work as long as the . _ : 
constellation distortion is not severe, so that the 
Rayleigh and Rice distributions are approximately .valid 
10 for their appropriate situations, 

The n-th moment of a probability distribution Px{x), --^ 
in the random variable .X, is here denoted,. _ 

/i;=Eixi, - (4) 

/i„ = E|(X-/iOl. <5) 

15 is the variance of the distribution. One measure that 

is useful for characterizing probability distributions 
is the ratio of the mean to the standard deviation, 

P.=/^./V^^. (6) 
Two other useful measures are the skewness, 

20 Py^l^Jif^.y. (7) 

and the kurtosis, 

Pa=/^.//^' (8) 



(see N. L. Johnson and S. Kotz, Distributions in 
25 Statistics, Vol. 1, Discrete Distributions , Houghton 
Mifflin Company, 1969, p. 18). The skewness is a 
measure of the asymmetry of the distribution. Positive 
values of p3 indicate that the distribution is "skewed to 
the right" so that the right tail is in a certain sense 
30 heavier than the left. Similarly, negative values 

indicate that the left tail is heavier than the right. 
Kurtosis is a measure of the peakedness of the 



wo 00/07302 



23 



PCTAJS99/16864 



distribution- For a normal or Gaussian distribution, p4 
= 3. For p4< 2, the distribution is less peaked than a 

Gaussian; f or p4 > 3, it is more peaked. 

Equations (6) to (8) are expressed in terms of 
5 central moments, except for fi{ since there is no Note 
that the shape factors py, pi, and can all be_ expressed 
in the form p^^/I^iptyY' , However, p2 is meaningless since ■ 

The central, moments can be escpressed in terms of 
10 the moments (4) by 



15 



20 



(9) 



For a Rayleigh distribution, as in (1), the 
moments (4) are given by 



/i; = (2cTT-r[i+^'j. (10) 

where r(k) is the Gamma function, with properties for 
k > 0, 

fkTfk) in general 

nk + l) = i . ■ ^ . (lla) 
t k! integers 

Inserting Equations (11) into (10) gives 



/i; = ^(r, {12a) 

/2;=2<T^, (12b) 

/i; = 3^<r\ (12G) 
^l\=^a*. (12d) 

25 Then, inserting Equatiohs (9) and (12) into Equations 
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(6) to (8) and simplifying gives the shape factors for 
the Rayleigh distribution as 

P.«^=^/7f^= 1-9131. - (IBa) 

24jt(n-3>) _ mK\ 
Pi^= (4-;r)^' =0-6311. (13b) 

P4b.^=-^^—^ =3.2451. (13c) 

For the Rice distribution in Equation (3), which is 
equivalent to a noncentral chi-square distribution with 
two degrees of freedom, the moments (4) are given by 

^; = (2a=fe-""T[f*l],F,(^^l;l;^). (14) 

10 where ,FI(o:,j3:X) is the confluent hypergeometric function, 
defined as 

^^r(a)r(^+k)k! P*"^^ • 

Comparison of Equations (10) and (14) shows that the 
moments of the Rice distribution are related to those of 
15 the Rayleigh distribution by 



Furthermore, defining 



' n , . cr 
2 2a' ) 



(15) 



2=— T, - (16) 

2a- 

we have that as z->0, e"^— >1 and ,FJ(cr,l;2) 1 for arbitrary 

20 a, so that ^^1^^^ fi'r^^^ and Equation (14) reduces to 

Equation (10) . Note that z is the interference-to-noise 
ratio (INRJ . 

Using published transformations of the confluent 
hypergeometric function, the moments of the Rice 
25 distribution may be expressed as 
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^i«ce = -^®"'[i Pi(-0-5;l;z) + 2z ,F,.(0.5;l;z)> = q(z)a, <17a) 
/^2«.„ = 2(1+ z)e-^F,(l;l;z)c7-= G(z)a-. 



/^BRice = -/f ® ' + Z) , F,(-0.5;1;Z) + (4Z^ + 8z- 1 ), F, (0.5; l;Z)p= C,{Z)a\ .- 

^4Ricc = 8e"',T, (3; 1;z)(t'* = 4(2" + 4Z+ 2)t'' = Q(Z)<t'*, (17d) ' 

where ,F,(-0.5;1;Z) , ,F,(0.5;1;Z) and ,F,(1;1;Z) are all tabulated 
functions (see M. Abramowitz and I. Stegun, Handbook of 
Mathematical Functions, Dover Publications, Inc., 1965, 
10 pp. 504-534) . Inserting these expressions into 

Equations (9) gives the central moments for the Rice 
distribution as 

n,(z) = -{q(z)-[qu)f]<y\ (i8a) 

= ■!C,(Z)-3Q(Z)q(Z) + 2[q(z)f }a\ (18b) 

)u,(z) = -|q(z)-4q(z)q(z)+6Q(z)iq(z)]'-3[q(z)]'y. (iso 

Then, inserting Equations (17a) and (18) into Equations 
(6) to (8) gives the probability distribution shape 
factors p„. The results are 

p,(2)=-=Sl] (19a) 

7C(Z)-[q(Z)]- 



20 



q(z)-3q(z)q(z)+2[q(z)f 
jc;(z)-[q(z)]^]' 



p,(2)= ^ AT ' , (19b) 



q(z)-4q(z)q(z)-f6q(z)rq(z)f -3[q(z)r 
|;(z)-[q(z)]-} 

Equations (19) are plotted in Figures 7 to 9. 
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The Rice distribution is known to approach the 
Rayleigh distribution as z->0, and to approach a Gaussian 
distribution for large z. Note in Figures 7 to 9 that 
the curves for />„ have the values given in Equations (13) 

5 for the Rayleigh distribution at z = 0 . For large values 
of z, the kurtosis, p^, approaches the limiting value of 
3, as for a Gaussian distribution: ; Similarly, the ' , 
skewness, p^, approaches 0, indicating that the 
distribution is becoming more symmetric, like a Gaussian 

10 distribution /"For small z/ 'p^ is large, indicating that 
the distribution is highly - asymmetric . 

It appears-.f rom Figures .1 to 9 that a reasonable 
test as to whether a CW interferer is present in the 
data, and corrupting the constellation, is to calculate 

15 the ratios p, , P3 , and , for the observed data averaged 
over all symbol clusters, and compare the values with 
Figures 7 to 9. If p, >110, P3<0.50, and p^<3.10, one can 
be reasonably certain that the radial probability 
distribution of points in the clusters is not Rayleigh, 

20 and that there is probably an interferer present. Based 
upon the curves in these figures, if a single parameter 
were to be used to make this decision, it should be the 
skewness, P3 , which exhibits a rapid and monotonic 
departure from the Rayleigh value at z = 0 . However, 

25 the kurtosis, P4 , also exibits a useful distinction 
between large values for small z (small or no 
inter ferers) and smaller values for large z (large 
inter ferers) , although the kurtosis function is not 
monotonic . 

30 To check the validity of this approach, a number of 

simulations were run for a range of input CNRs and a 
range of interferer amplitudes, where the amplitude a 
was normalized to the noise margin d, and the normalized 
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interferer amplitude was varied over the range 0.000001 
< a/d < 0.5.-. Although most runs were done with- an 
amplitude^^^spread of 0 on the interferer (no .^amplitude 
modulation), some runs were also done with a^uniformly 

5 distributed amplitude spread of 0.1 to represent a 
modest amplitude distribution. Similar results were 
obtained with and without the amplitude modulation. The 
results of these simulations showed a strong correlation 
between the 'cluster skewness P.r and the calculated 

10 carrier-to-interference ratio (CIR) , and a' somewhat 

weaker correlation between the cluster kurtosis and 
the CIR.-' The distribution shape parameter Pi showed a 
much poorer correlation with the CIR. This prediction 
that the cluster skewness is the best single predictor 

15 of the presence of an interferer was subsequently 

confirmed by tests on a developmental VDSL modem that 
had the capability to make the various measurements and 
tests. Consequently, although both the skewness and 
kurtosis are shown in Step 54 of the procedure in Figure 
2, the skewness is the preferred parameter if just one 
is to be calculated. 

To calculate the skewness in Step 54, the symbol 
coordinates each sample symbol in each cell m are 

adjusted by subtracting the sample mean coordinates 

25 and y^^ to form the coordinate differences from the 
cluster center, ^k"^ ymk~ym- Then the radial 

distance of each sample from its cluster center, 

^mk= V (^mk" ^ + (Vrnk " Vm)^ • 
is calculated. If there is no distortion of the 
30 constellation, so that the expected values of the sample 
mean coordinates and coincide with the ideal 
locations of the cluster centers, then the sample means 
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in Equation 20 may be replaced by the ideal cluster 
coordinates. This also eliminates statistical^-'^ 
fluctuations in - the -calculated cluster-mean coordinates. 
Next, within' Step 54, 'the sample- moments 

5 .. _sJi = — XZ^mk. n=l,.,.,4, - 

are calculated, where is the number of data symbol 
samples in the bin for symbol cell m, _ accumulated 
previously in Step 48, and N is the total number of data 
symbol samples in all bins. Analogous to Equations (9), 
10 the second and third central moments are then calculated 
by ' ' • ^" 

&> = s;-s,'-, 
Sj = s^ -3s^s[ + 2Si'^ 
The cluster skewness is then calculated by 

15 p3 = ^/ihf' 

Similarly, to calculate the kurtosis in Step 54, the 
fourth central moment is calculated by 

s^=s;-4s^s,' + 6s;s,'^-3s[\ 
and the cluster kurtosis is then calculated by 

20 P4 = Si/s2. 

To determine whether an interferer is present, the 
calculated skewness and/or kurtosis is compared, in 
Step 56, with a threshold level determined from Figure 8 

25 or 9, respectively. As seen from Figure 8, a skewness 
less than about 0.50 indicates that the cluster radial 
distribution is Ricean distributed, indicating that an 
interferer is present, while a skewness greater than 
this value indicates a Rayleigh distribution, indicating 

30 that no interferer is present. Similarly, from Figure 9 
it is seen that a kurtosis less than about 3.15 
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indicates that the cluster radial distribution is Ricean 
distributed, indicating that an interferer ,.is_ Present, 
while a kurtosis greater than this value^.iiidicates a 
Rayleigh distribution, indicating, that _„,no _inter-ferer is 
5 present. In experiments on the developmental VDSL 

modem, a skewness threshold, of - 0 . 3 to 0 . 5 and a kurtosis 
threshold of 3.12 to 3.15 have been found to be optimum. 
However, the lower the CNR of the input signal, i.e., 
the higher the noise ~level,^ the stronger the i'nterferer 
10 must generally be to detect it above the noise. Of' * 
course, the smaller the interference-to-noise ratio 
(INK) is, the less" 'pro~poftionate effect "the'^interferer 
has on the overall SER. When the cluster skewness is 
plotted as a function of both the CNR and the normalized 
15 interferer amplitude, as in Figure 10 (which was 

generated for a 4-QAM signal), it is found that, except 
for some anomalous behavior in the vicinity of CNRs from 
15 to 20 dB (which. is reproducible, but of undetermined 
cause) , a given skewness threshold level can only 
20 detect larger amplitude interferers as the CNR 

decreases. Similarly jagged and anomalous results have 
been found when using the cluster kurtosis as the 
interference detection measure, as shown in Figure 11. 
With either approach, neither the skewness nor the 
25 kurtosis can be used to determine the interferer 's 

amplitude. But either can be used to determine whether 
an interferer with amplitude significantly (more so for 
kurtosis than for skewness) higher than that 
corresponding to that of the test criterion is present 
30 (except at low input CNRs and high interferer amplitudes 
where the contours of skewness or kurtosis versus CNR 
are essentially vertical) . Within these limits, for any 
CNR, a skewness and/or kurtosis threshold can be 
selected to determine whether an interferer with 
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amplitude greater than a specified level is present. 
For high noise and strong interf erers.,: the combined 
effect is to expand the sizer.of the symbol clusters-so .r- 
much that they overlap, - at which .point /determination of 

5 CNR, SER, and cluster skewness and kurtosis from . ■ ~ 
constellation analysis is no longer possible. Also, _ 
when there is much noise (low CNR) • the resulting - • 
statistical fluctuations, cause some -uncertainty . in the _ 
calculated skewness and -kurtosis,. reducing the - Cfi 

10 significance of .the. threshold -decision .\ucNevertheless , 

at moderate and high CNRs, the method .works quite well,. ^ 
and has been found valuable • in detecting weak . t/"- 
interferers. - - ^ 

Note from Figure 7 that the quantity 'p] = ^/^Js^ has a 

15 unique value for any value of the quantity z given in 
Equation (16), which is the INR. Therefore, if it has 
been concluded that an interf erer is present, the value 
of Pi may be calculated, and the corresponding value of z 
may be looked up in a tabulated version of Figure 7 to 

20 obtain an estimate of the INR for diagnostic purposes. 

When the skewness and/or kurtosis test 56 has found 
the presence of an interf erer, the computer 7 6 sends a 
notification message via reverse control channel 58 to 
transmitter controller 60. Transmitter controller 60 

25 then interrupts the input bit stream 20 and turns off 
the power amplifier 28, so that the only input to the 
receiver 18 is the external interference 32. Meanwhile, 
the computer 76 also sends a notification message to 
receiver controller 62, which sets the carrier frequency 

30 arid symbol rate synchronizers to flywheel, opens both 
poles of switch 64, and closes switch 70, so that the 
received external interference 32 goes (through an 
amplifier, not shown) ^to the analog- to-digital 
conversion 78 and spectrum analysis 74. The receiver 
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controller 62 also sends an acknowledgment to the 
transmitter controller 60 that the receiiver 18 is - 
temporarily ceasing to process -received digital data 
while the computer- is performing the spectriim analysis. 

5 The transmitter controller 60 does not turn off -the 
transmitter 16 until this aclcnowledgment is received. 
Simultaneously, the computer 76 switches its operation 
entirely over to -spectrum analysis- 74.-"The results of 
the spectrum analysis are sent via reverse control - 

10 channel 58 to transmitter controller '60 which analyzes 
the received spectrum to find a carrier frequency at 
which the desired'data rate may be transferred without 
encountering any interference. If no pass band of 
sufficient width is available, the transmitter 

15 controller 60 will switch transmission to a higher order 
constellation with a lower symbol rate, providing that 
such a combination is available that can provide the 
desired bit rate, and provided that the noise margin in 
the new constellation is adequate to provide the 

20 required BER. If these conditions cannot be met, the 
transmitter controller 60 will attempt to reduce the 
transmitted bit rate with the same or larger 
constellation in order to reduce the transmitted 
bandwidth, providing that the information to be 

25 transmitted will not be impaired by a reduced bit rate, 
as would some types of real-time transmissions. 

After selecting a new carrier frequency, 
constellation size, and/or symbol rate, the transmitter 
controller 60 inserts a control message into the input 

30 bit stream 20 to transmit over the forward transmission 
channel 3 0 to the receiver 18. The receiver controller 
62 extracts this message out of the output bit stream 4 6 
and then sets the receiver 18 to receive and demodulate 
a signal with these neV parameters. The receiver 
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controller 62 also sends an acknowledgment of the 
receipt of -this control message over -.the reverse . control 
channel 58- to the .transmitter controller .60, -Only after - 
receiving .this acknowledgment does -the transmitter - 

5 controller 50 actually change the transmitter parameters 
to the new -.values .... 

At this point, receiver controller 62 also opens ^ 
switch 70, closes both poles of switch 64, and directs 
computer_76 to switch back to accumulating statistics * on 

10 the new received symbol streams calculating the CNR and 
BER; testing ^ the. BER; and, if necessary; calculating and 
testing the skewness and/or kurtosis on the new received 
constellation. - - - = 
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CIAIMS: " "'^^ ■ 

1. In a communication apparatus"/ a method of detecting- 
5 whether an interference" signal is being received 

with a received signal of interest (SOI), said 
method comprising steps of: 

' (a) accumulating statistical inf ormation— - 

associated with said SOI;* '■ - ~ ' 

10 (b)' determining whether 'ah intrusive - 

interference" signal is present with-said 
SOI; '~ ' ' ' - - - 

(c) analyzing, when an intrusive interference 
signal is present, said intrusive 

15 interference signal to determine the 

frequency thereof; 

(d) identifying a spectral region free of 
interferers; and 

(e) adjusting said communication apparatus to 
20 operate in said spectral region 

identified in step (d) . 

2. The method according to claim 1 wherein said steps 
(a) through (d) are non-intrusive to a continued 
25 reception of said SOI in said communication 

apparatus . 



3. The method according to claim 2 wherein said step 
30 (b) further comprises steps of: 

(f) determining a root mean square (RMS) 
cluster size of said intrusive 
interf ei^ence signal; 

(g) determining a reciprocal constellation 
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10 



closure of said intrusive interference 
„ signal; and 
(h) determining a carrier-to-noise ratio 
. (CNR) of said intrusive interference - 
signal- 



The method^. of claim 3 wherein said step (b) further 
comprises determining a syrrfDol error rate of said 
intrusive interference signal. 

The method of claim 3 wherein said step (b) further- 
comprises determining a bit error rate of said 
intrusive interference signal. 



15 6. The method of claim 3 wherein said step (b) further 
comprises determining an average cluster skewness 
associated with said intrusive interference signal. 

7 . The method of claim 3 wherein said step (b) further 
20 comprises determining an average cluster kurtosis 

associated with said intrusive interference signal. 

8. A transceiver for transmitting and receiving data 
signals, said transceiver being capable of non- 
25 intrusively detecting an interference signal to said 
data signals, said transceiver comprising: 

storage means for accumulating statistical 
information associated with said data signals; 

interference determining meeins for determining 
30 whether an unacceptable interference signal is present 
with said data signals; 

analysis means for spectrum analyzing said 
unacceptable interference signal and determining a 
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frequency thereof and further locating a spectral region 
free of inter ferers; and — - i c 

conveyance means for communicating to a - 
transmitter said spectral region -fre.e of interferers- for 
5 changing operation of the transceiver and transmitter to 
said spectral region free of interferers. 

9. The transceiver ^of claim 8 .wherein said- 
interference signal is a continuous wave signal. 

10 0 

10. The, transceiver of claim.- 9 ^.wherein, said transceiver 
uses quadrature amplitude modulation. . : 

11. The transcevier of claim 9 wherein said transceiver 
15 uses carrierless AM/PM data communication methods. 

12. The transceiver according to claim 9 wherein said 
storage means acciimulates both presliced quadrature 
amplitude coordinates and sliced symbol cell 

20 identification data. 

13 . The transceiver according to claim 12 wherein said 
analysis means processes a sequence of demodulated 
symbols to calculate the statistics of in-phase and 
25 quadrature mean coordinates and in-phase and quadrature 
sample variances for a subset of symbols in each 
decision cell from said sequence of demodulated symbols, 
said sliced data being used to assign said presliced 
data to a proper symbol cell. 

30 

14. The transceiver according to claim 13 wherein said 
analysis means further determines from said statistics a 
root mean square (RMS) cluster size for all the symbol' 
clusters in a constellation, a constellation closure, a 
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resulting carrier-to-noise-ratio, a resulting symbol 
error rate,, and a resulting _bit error rate (BER) . 

15. The- transceiver according to Claim 14 -wherein said 
5 analysis means further tests the magnitude of said bit 

error rate to determine whether the quality^ of -said 
demodulated signal is above an acceptable threshold and 
if not,:- then directing the commencement of. further" 
processing to determine whether .the. degradation of the 
10 quality of said demodulated signal is due to the 

presence of said continuous wave interfering .signal . 

16. The Transceiver according to claim 14 wherein said 
analysis means determines an average cluster skewness 

15 for the symbol clusters in said constellation and tests 
the magnitude of said skewness to decide whether an 
interferer is most likely present, which would be the 
cause of said degradation of said BER. 

20 17 . The transceiver according to claim 14 wherein said 
analysis means determines an average cluster kurtosis 
for the symbol clusters in said constellation and 
testing the magnitude of said kurtosis to decide whether 
an interferer is most likely present, which would be the 

25 cause of said degradation of said BER. 

18. The transceiver according to claim 14, wherein 
said conveyance means transmits the results of said test 
of skewness back to said transmitter over a reverse 
30 channel, said transmitter temporarily turning off 
transmission while said receiver receives only the 
external interference signals in the RF environment, and 
spectrum analyzing said^ environment to determine the 
frequencies of said interference signals. 
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19. The- transceiver according- Jtq- claim -14 wherein said 
conveyance means transmits the results of said test of 
kurtosis back.. to^said rtransmitter- oyer ~a, reverse 

5 channel, said transmitter temporarily turning off 
transmission while said receiver .receives only the 
external interferencer signals in the RF -environment, and 
spectrum analyzing said environment - to determine the 
frequencies, of - said interference signals. 

20. The- transceiver according to claim 14. wherein , said 
conveyance means transmits the results of said spectrum 
analysis back-_to said transmitter over a reverse 
channel, said transmitter determining new transmission 

15 parameters to avoid said interference signals, and said 
transmitter transmitting said new transmission 
parameters to said transceiver. 

21. The transceiver according to claim 20 wherein said 
20 transciever and said transmitter reconfigure to transmit 

and receive the desired data signals with said new 
transmission parameters and then resuming data 
transmission using said new transmission parameters, 

25 

22. The transceiver according to claim 20 wherein said 
new transmission parameters include a new carrier 
frequency, QAM modulation order, and data rate. 

30 23." The transceiver according to claim 9 wherein said, 
transceiver uses quadrature part^ial response modulation. 

24. A modem capable of non- intrusively detecting an 
unwanted interferer during a receiving of a data signal, 
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comprising: 

a transmitter for ' transmitting information to a 
communications apparatus ; ^ - 

a receiver- for receiving said data- signal from said 

5 communication' apparatus; - - . - 

a memory coupled to said receiver-- for storing : . ^ 
statistical information '-associated with' said data* . : 
signal; and- ^ - - ^- ^ „ 

a processor for measuring ^-constellation' statistics 

10 from said statistical information, determining a 

carrier-to-noise" ratio and bit error--rate of said data . - 
signal, determining whether -a continuous- wave or - 
narrowband interferer is present, determining whether 
said bit error rate exceeds an acceptable limit by 

15 testing whether a symbol cluster changes from a Rayleigh 
distribution to a Ricean distribution, and when said bit 
error rate exceeds said acceptable limit, interrupting 
said received data signal to spectrum analyze a 
transmission band for an interference free region of 

20 said transmission band and communicating said 
interference free region to said communications 
apparatus for adjustment. 

25 25. The modem according to claim 24 wherein said 

testing comprises calculating a skewness associated with 
said symbol cluster. 

26. The modem according to claim 24 wherein said 

30 testing comprises calculating a kurtosis associated with 
said symbol cluster. 

27. The modem according to claim 23 further comprising 
an analog-to-digital converter coupled for converting . 
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said data signal into a digital signal. 

28. The modem according to claim 25 further comprising 
a spectrum fast-Fourier-transform-analysis^ means for- • 
5 spectrum analyzing said analog signal. 
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